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ABSTRACT 



We present the gravitational wave analysis from rotating (model sl5g) and nearly non-rotating (model sl5h) 3D MHD core collapse 
supernova simulations at bounce and the first couple of ten milliseconds afterwards. The simulations are launched from ISMq pro- 
genitor models stemming from stellar evolution calculations. Gravity is implemented by a spherically symmetric effective general 
relativistic potential. The input physics uses the Lattimer-Swesty equation of state for hot, dense matter and a neutrino parametrisa- 
tion scheme that is accurate until the first few ms after bounce. The 3D simulations allow us to study features already known from 
2D simulations as well as nonaxisymmetric effects. In agreement with recent results we find only type I gravitational wave signals 
at core bounce. In the later stage of the simulations, one of our models (sl5g) shows nonaxisymmetric gravitational wave emission 
caused by a low T/\W\ dynamical instability, while the other model radiates gravitational waves due to a convective instability in the 
protoneutron star. The total energy released in gravitational waves within the considered time intervals is 1.52 x 10"'' Mq (sl5g) and 
4.72 X IO^'^Mq (sl5h). Both core collapse simulations indicate that corresponding events in our Galaxy would be detectable either 
by the LIGO or Advanced LIGO detector 
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1. Introduction 

] Gravitational wave astronomy could well make soon its first ob- 
servations as the running earth-based facilities such as LIGCQ, 
VIRGCB GEO6OC0, TAMAQ and AIGCS are up to reach the 
required high sensi tivities. For a review o n gravitational wave 
detection see e. s. (iHough & Rowanll2007i) . The successful ob- 
servation of gravitational waves (GW) would be a major break- 
through: it would open a new window to the universe, allow- 
ing u s to observe electromagnetically hidden regions (Thornd 
I Il995h . One of the most promising gravitational wave sources is 
, the stellar core collapse leading to a supernova explosion. The 
observation of both, GW and the neutrino signal, from a galactic 
supernova would reveal hitherto hidden details about the explo- 
i sion scenario and impose constraints on the nuclear and weak 
' interaction physics under conditions that can not be obtained in 
terrestrial experiments. 

As a core of a star of M > 8M0 reaches the end of its stellar 
evolution, it becomes gravitationally unstable as soon as ther- 
monuclear burning produces a significant amount of iron group 
nuclei. Next to photo-disintegration, electron captures on free 
protons and nuclei reduce the mostly electron supported pres- 
sure, and the core starts to collapse eventually. The collapse con- 
tinues until nuclear densities of ~ 2 x 10'''g/cm^ are reached, 
depending on the equation of state (EoS). As soon as we enter 
this density regime, the core overshoots its equilibrium position 
and bounces back. A sound wave immediately forms and steep- 
ens into a shock wave that propagates outwards. However, the 
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shock stalls within ~ 5 ms after core bounce due to the large 
energy loss caused by the dissociation of nuclei into free nucle- 
ons at a cost of ~ 8.8 MeV per nucleon and by neutrino emission 
connected to copious electron captures on the emerging free pro- 
tons. It continues to propagate outwards to radii around 100-200 
km as standing accretion shock. Hereafter, i.e. some ms after the 
prompt explosion mechanism failed, a delayed explosion mech- 
anism by neutrino heating is thought to occur, e.g. Bethe ( 19901) . 

The idea of reviving the stalled shock again via neutrino 
reactions behind and ahead of the shock ha s long been in- 
vestigated as possible e xplosion mechanism jBethe & WilsonI 
ll985VlJankaetani2001h . In the dissociated matter behind the 
shock, electron-flavour neutrino captures are the dominant heat- 
ing reactions, while in the accreting matter ahead of the shock, 
electron-neutrino absorption and neutrino-nucleon scattering are 
the domi nant interactions that pr eheat the infalling unshocked 
material (iBruenn & Haxtonll99lb . In addition, neutrino-electron 
scattering and pair annihilation may contribute to the heating as 
well. However, the effects of pre-heating and neutrino-electron 
scattering on the shock revival are insignificantly small com- 
pared to the effect from the capture of neutrinos on free nucleons 
in the heating region below the shock. Equally as important as 
the neutrino heating is the neutrino cooling of matter that settles 
on the protoneutron star (PNS) ( Janka 2QQ%. 

As a core collapse superno va is likely to show aspherical 
features also close to its center (' Leonard et al.ll2066h . where the 
matter assumes a high density, a tiny fraction of the released 
binding energy can be emitted via gravitational radiation. The 
following features have been suggested as possible causes of 
the asymmetries in the energy-matter distribution necessary to 
emit GW: the rotational stellar core collapse, convection in the 
high-density protoneutron star, nonaxisymmetric rotational in- 
stabilities, fluid instabilities in the lower-density hot mantle sur- 
rounding it (possibly triggering neutron star oscillation modes). 
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and an anisotropic neutrino emission (fo r recent reviews see e.g. 
(iKotake et al.ll2006HFrver & Newl l2006')). 

The understanding of GW emission from differentially ro- 
tating core collapse has evolved with time due to the improv- 
ing input physics use d in sim ulations. In 2D axisymmetric com- 
putations by MiilleB (Il982l) iron cores, a Newtonian hydrody- 
namics code and a tabulated finite temperature EoS were used. 
The results allowed to recognise the link between rotation and 
the efficiency of GW emission. After having taken into ac- 
count more micro-physics in Newtonian gravity, such as electron 
capture on protons a nd a simplified neutrino transport scheme 
jMoenchmever et al] IT99 1 ). one was able to distinguish two 
dif ferent wave form charact eristics, later known as type I & 
II. IZwerger & Miilleii (Il997h performed an extensive parame- 
ter study of a wide variety of models, using progenitors in ro- 
tational equilibrium, Newtonian gravity, a polytropic EoS, but 
neglecting electron capture and neutrino physics. The results 
lead to the introduction of type I, II & III wave forms and 
their quantitative distinction in dependence of the stiffness of 
the EoS and rotation. A type I waveform is characterised by a 
large amplitude peak at core bounce and subsequent damping 
ring-down oscillations. It appears if the influence of the initial 
angular momentum is small. In this case the core bounce oc- 
curs by the stiffening of the equation of state and is pressure- 
dominated. A type II signal occurs if the core bounce around 
nuclear density is dominated by strong centrifugal forces; it 
has several distinct peaks caused by multiple 'centrifugal' core 
bounces with following coherent re-expansion phases of the in- 
ner core. The type III signature is defined by a 'large' positive 
peak at bounce followed by some smaller oscillations with very 
short periods. It appears in the case of fast, pressure-dominated 
core collapse due to very efficient electron capture. Continued 
improvement of the input physics by using realistic equations 
of state or neutrino physics, by adding magnetic fields or in- 
corporating GR effects, and by extending the dimensionality 
from 2D to 3D, permitted to deepen and underline the quali- 
tative and quantitative understanding of GW forms frorn rota- 
tional core collapse around bounce time (Janka & Muller"l996 ; 
Rampp et al. 1998; Frver et al. 2002; Dimmelmeier et al. 2002; 
Kotake et al. 2003, 2004; M uUer et al j2004i; IDimmelmei er et all 
2005i, ,2007; .Ott et al...2007|). In the most recent simulations by 
Ott et al.1 (l2007l) and iDimmelmeier et all (|2007|) it was pointed 
out that realistic input physics (GR, a micro-physical EoS, ap- 
proximative description of deleptonisation) might lead solely to 
type I wave forms. 

The investigation of GW emission resulting from convec- 
tively driven small-scale aspherities inside the PNS, the 'hot- 
bubble' and anisotropic neutrino emission are still hampered by 
the requirement of computationally expensive neutrino transport 
in the postbounce evolution. Nevertheless, simulations that in- 
clude accurate neutrino transport in axisymmetric models were 
explored with respect to GW in Miilleret al. (2004) based on 
state-of-the-art progenitor models, GR-corrections, and sophis- 
ticated equations of state. The postbounce phase is believed to 
provide for hundreds of milliseconds, or even longer than 1 s, an 
interestingly large signal with a frequency distribution between 
~ 300 - 1200 Hz in case of convection, while the dominant fre- 
quencies lasting from anisotropic neutrino emission lie at about 
Si 10 Hz. Nonaxisymmetric dynamics has rec ently been inves - 
tigated with neutrino physics approximations (lOtt et al.ll2007l) . 
In these computations it was possible to show without the ad- 
dition of any seed perturbation that differentially rotating pro- 
toneutron stars are likely to become dynamically unstable at low 
j3 - r/|W|-values (ratio of rotational to gravitational energy). 



leading to s t rong ii a rrow-band GW e m ission (ISaiio et al. (20031); 
Watts et al.' ('200f) jQttetiD (l2005h ; ISaiio & Yoshidal (l2006h ; 



Ou&Tohline (20061)). 

Another mechanism for gravitation al wave emission in core 
collapse supernovae was proposed in lOtt et al.l (l2006l) . It was 
pointed out that PNS core g-mode oscillations might be a major 
source of gravitational radiation. The release of GW in their sim- 
ulations is highly correlated with the fundamental core g-mode, 
leading to a spectrum that peaks at twice the frequency of the 
aforementioned g-mode. 

In this article we present the gravitational wave analysis 
from two 3D magneto-hydrodynamics (MHD) core collapse 
simulations. The simulations are based on a nuclear equation 
of state. The deleptonisation by electron capture and neutrino 
emission are accurately parametrised up to about 5 ms after 
bounce. Important general relativistic corrections to the gravi- 
tational potential are taken into account as well. The simulations 
are launched from progenitor models at the end of stellar evo- 
lution calculations (IWooslev & Weaveill 1995b . One of the initial 
models (sl5g) rotates moderately fast at the onset of the collapse 
with 27r rad/s for the central angular frequency, whil e the second 
model (sl5h) rotates more slowly as estimated in iHeger et al.l 
(12005^. We perform contiguous simulations from collapse into 
the postbounce phase and study three-dimensional effects on the 
emission of gravitational waves. The GW-producing asymme- 
tries in three-dimensional simulations could be significantly dif- 
ferent from the asymmetries present in axisymmetric models. As 
only two other group have been performing three-dimensional 
supernova models with a nuclear equation of state and neutrino 
physics approximations for the prediction of gr avitationa l wave 
signals (Fr yer et al. 20041; iDinimelmeier et al.l 2007; O tt et al.l 
2007!) our simulations provide an independent confirmation of 
the latter results. Rather than augmenting the database with nu- 
merous GW pattern templates, we aim to underline the impor- 
tance of the neutrino physics in supernova GW signal predic- 
tions and to encourage the search for GW signals from Galactic 
core collapse supernova events. 

The paper is organised as follows: In Sect. II we define the 
numerical methods and input physics used in our supernova 
models. We also describe the gravitational wave extraction for- 
malism applied for the analysis of our three-dimensional data 
sets. Section III is devoted to the resulting gravitational waves, 
their properties and the possible detectability of signals. In Sect. 
IV we discuss the uncertainties in the numerical models and try 
to separate the robust features from details that are subject to 
change in future improved models. A conclusion is given in Sect. 
V. 



2. Method 

2.1. Description of tine liydrodynamical models s15g and 
s15h 

Only few recent multi-dimensional collapse simulations made 

the effort to include neutrino physics. These schemes are either 

i 1 1 1 

very computationally expensive (Buras et al. 2003; Dessa rt et al.l 

2006) or rely on simplifications of the neutrino trans port an d 
its micro-physics (iKotake et all 120031: iFrver & Warrenll2004al) . 
The complete Boltzmann neutrin o transport equatio n can only 
be solved in spherical symmetry (( lMezzacappal2005h and refer- 
ences therein). Hence, three-dimensional hydrodynamical mod- 
els must rely on neutrino transport approximations. 

A simple and computationally efficient parametrisation of 
the deleptonisation in the collapse phase can be derived from 
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Fig. 1. Comparison of the almost non-rotating 3D model sl5h 
with the spherically symmetric model G15, which is based on 
general r elativistic three-flavour Boltzmann neutrino transport 
jLiebendorfer et al. 2005) . From the upper left to the lower right 
we compare as a function of enclosed mass: the density-, the 
velocity-, the Yg-, and the entropy profiles. The solid lines show 
the results of model si 5h and the dotted lines the results of model 
G15. The thin lines represent a time instance at 5 ms before 
bounce and the thick lines represent a time instance at 5 ms after 
bounce. Excellent agreement is found in all four quantities — 
with one exception: The parametrised neutrino leakage scheme 
cannot model the neutrino burst. The neutrino burst causes a 
prominent y(,-dip and additional cooling in the G15 data. 



a tabulation or fit of the electron fraction as a function of 
density. S pherically symmetric mod els with Boltzmann neutrino 
transport jLiebendorfer et alj|2005l) show that the Yg in different 
layers in the homologously collapsing core follow a quite sim- 
ilar deleptonisation trajectory, i.e. reach similar Y^ values as a 
function of density. As changes in Y,, can only be due to elec- 
tron captures, it is possible to deduce corresponding changes 
in entropy an d to calculate the n eutrino stress from the emit- 
ted neutrinos ( lLiebend6rfeiJl2005 '). Figure [T] shows a compari- 
son of the 3D parametrised run sl5h, launched from an almost 
non-rotational progenitor model, with a spherically symmetric 
general relativistic model G15 th at is based on general rela tivis- 
tic Boltzmann neutrino transport dLiebendorfer et all2005h . The 
ID results of core collapse are accurately reproduced by the 
3D run. The parametrised neutrino physics presents a signifi- 
cant improvement with respect to adiabatic simulations and may 
even rival with neutrino transport schemes that neglect neutrino- 
electron scattering. However, the accuracy breaks down with the 
launch of the neutrino burst at a few milliseconds after bounce. 
With the currently implemented scheme, accretion flows in the 
postbounce phase deleptonise only down to Y^ ~ 0.3 instead of 
Ye ~ 0.15. Neutrino heating is neglected altogether. This turns 
the quantitative model of the collapse phase into a qualitative 
model of the postbounce phase. 

The hydrodynamics of the 3D simul ations is based on a sim- 
ple and fast cosmological MHD code jPen et alj|2003h which 
has been parallelised, improved and adapted to the supernova 
context. A realistic equation of state (Lattimer & Swestv 1991) 
is used and the monopole term of the gravitational potential is 
implemented by a spherically symmetric mass integration which 



includes general relativistic corrections dMarek et al.ll2006h . The 
3D computational domain consists of a central cube of 600^ km-' 
volume, treated in equidistant Cartesian coordinates with a res- 
olution of 1 km. The 3D code has buffer zones that require at 
each time step a prescription of the conditions in the immediate 
neighborhood of the 3D computational domain. In order to ob- 
tain these conditions, we embed the 3D computational domain 
in a larger spherically symmetrical computational domain that is 
treated by a spherically sy mmetric hydrodynamics code 'Agile' 
jLiebendorfer et all2002 i). After each time step, the buffer zones 
at the boundary of the 3D hydrodynamics code are filled with the 
current conditions of the spherically symmetric solution. As long 
as the infall velocity at the boundary is subsonic, we spherically 
average the conditions in the 3D code and feed them back to the 
mass shells of Agile that are enclosed by the 3D domain. Once 
the infall velocities become supersonic, this step can be omitted 
because no hydrodynamic signal can leave the 3D domain any- 
more before the shock passes the boundary after the onset of an 
explosion. 

For both of our mode ls we use a 15 M© progenitor star of 
jWooslev & Weaveill 19951) as initial model. As the rotation rates 
of inner stellar cores are not very well-known, we assign angular 
momentum according to a simple parametrisation with a shellu- 
lar quadratic cutoff at 500 km radius. The angular momentum is 
assumed to be conserved until the infalling layers enter the 3D 
computational domain. For one model, which we name sl5h, we 
set an initial angular velocity of O = 0.3 rad/s as obtained in a 
stellar evolution model t hat takes account o f the effects of rota- 
tion and magnetic fields jHeger et al.ll2005i) . Because this value 
represents a rather slow rotation, we performed an alternative 
model, sl5g, with an initial angular velocity of Q = 27r rad/s. 
This is still a slow rotation rate compared to values assumed 
in many parametrised studies of the GW signal, but it allows 
for a clear distinction from model sl5h. The initial - T/\W\ 
comports 0.26% (sl5g) and 0.00059% (sl5h), respectively. Both 
models collapse in a similar manner to spherically symmetric 
simulations. After bounce, the declining strength of the bounce- 
shock results in a negative entropy gradient and the protoneutron 
star becomes convectively unstable on a short time scale. The 
simulation sl5g was carried out until about 70 ms postbounce, 
whereas sl5h run ~ 100 ms postbounce. 

We h ave set the magnetic fields to the values suggested in 
i Heger et al. (2005). The advection of the magnetic field into the 
boundary of the 3D computational domain is technically very 
delicate because the field lines assigned to the buffer zones have 
to ensure a vanishing divergence at the transition to the slightly 
evolved magnetic field on the interior of the boundary. We found 
a quite pragmatic solution where the field is first assigned ac- 
cording to above analytical setup of a divergence-free initial 
field. In a second step, the divergence at the boundary is anal- 
ysed and corrected such that the transition is divergence-free. In 
any case, due to its weakness, the magnetic field does not influ- 
ence the hydrodynamical evolution noticeably during the early 
postbounce evolution. 



2.2. Extracting gravitational radiation in fhree dimensions 

Throughout this article we work in cgs-units and use the follow- 
ing values for the speed of light, c= 2.997 ■ 10'" cm s"', the 
gravitational constant, G= 6.672 ■ lO^'^cm^g^'s"^, and the par- 
sec, 1 pc = 3.086 -10'** cm. We do not assume any symmetry. 
In this case, the gravitational wave field hj^ can be decomposed 
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into two orthogonal polarisations with amplitud es A+ and Ax, 
see e.g. ( iMisner et al.lll973l:lFinn & Evanslll990h : 



and for 6 



hJ/{X,t)^-iA^e^+Ay,e^), 



(1) 



where R is the distance to the source. In spherical coordinates, 
the unit polarisation tensors are given by: 



and 

ex = - 



(2) 



(3) 



In the slow-motion hmit dMisner et alj|1973t iFinn & Evans! 

1 19901) the amplitudes A+ and Ax are given by linear combina- 
tions of the second time derivative of the transverse traceless 
mass quadrupole tensor: 



- tgg - f^. 

Ax — 2fflA. 



(4) 
(5) 



In the Cartesian orthonormal basis, the quadrupole tensor is 
given by 



J dVp ^XiXj - ^ 



3^ 



(6) 



Note that its time derivatives are evaluated at a retarded time. 
However, there exist shortcomings related to numerical high- 
frequency noise and the momentum-arm which make the per- 
formance of a direct evaluation of Eq. ^ poor, as discussed 

J " 

in (IFinn & Ev ans 1990). Therefore, we use an alternative post- 
Newtonian expression from Blanc he let alj (|T990), where the 
second order time derivatives of the quadrupole moment are 
transformed into hydrodynamical variables which are known 
from the core collapse simulation. 



j.TT _ ^ 



J dVp (iviVj - Xidj^ - xjdi<i>) , 



(7) 



where O is the gravitational potential. This expression allows 
the evaluation of the wave field from data at only one time in- 
stance; moreover, the integral has a compact support. Note, how- 
ever, that this formula is not gauge invariant and only valid in 
the Newtonian slow-motion limit. Nevertheless, it was shown 
that this approximation is sufficiently accurate when compared 
to other approaches (Shibata & Sekiguchi 2003). 

The polarisation modes can explicitly be obtained from a co- 
ordinate transformation, for example for the 00-component: 



jj dx' dx-i 



(8) 



This leads to the following non-vanishing components 

(lOoharaet al.lll997h : 

tgg - (f[J cos^ -I- sin^ + It^J sin (p cos (p) cos^ 9 



I sin cos 6 



+ sin^ 6 - 2{tlJ cos + fj/ sin ^ , 

- f sin^ 4> + t^J cos^ (p - 2?!/ sin ^ cos 

- (t^J^ - f ) cos 6 sin (pcos(l> + t^^ cos 6i(cos^ (p 



sm 



^ d)) + fj sin 6 sin ( 



• t^.I sin cos ( 



We evaluate below for simplicity only the gravitational wave am- 
phtudes for 9 - (p -Q (denoted below with the subscript /) 



A + l ?v-v tv- 



ix/ 



= 2f 



xy 



(9) 
(10) 



= (denoted as II) 



A +11 
Axil 



(11) 

(12) 



We notice that A+// corresponds to the non-vanishing quadrupole 
amplitude Af^ of axisymr netric rnodels in the multipole expan- 
sion of the radiation field (lThornelll980h . 

The energy carried away by gravitational radiation can be 
calculated by the following expression: 



5G J dt 

" TsgJ 



Hij - lSijIii)]^dt 

dt[h^ +Iyy +1^ 



(13) 



^.v.vA'v ^xx^zz A'v^z; 



where /,y - 

Whenever the question is raised whether the emitted signal 
of a source would possibly be detectable by earth-based facil- 
ities, such as e. g. LIGO, one is in particular interested in the 
detector dependent characteristical frequenc ies fc and ampli- 
tudes he (calculated according to Eq. (31) in (iHawking & Israeli 
1989)) and the 'signal-to-noise ratios' (SNR). We estimate the 
signal-to-noise ratios for optimal filte ring searches as follows 
(iFlanagan & Hugheslll998HMiiller et al..,2004.) : 



(-) 

\n). 



optimal/ titer 



f 



\h(y) 
Sh(v)' 



-dv 



\h{v)\v 



1/2 



d log V, 



(14) 
(15) 



where h{v) is the Fourier transform of the gravitational wave am- 
plitude, h{v) - h+{v) + hxiv), and S h(v)[l / Hz] is the power spec- 
tral density of strain noise in the detector and hrmAHz'^^^] = 
-\JSh(v) stands for the rms noise level of the detector at a given 
frequency v. The Fourier spectra were normalised by using 
Parseval's Theorem. 

Below, we determined the signal-to-noise ratios for an op- 
timal orientation of detector and source. As detector sensitiv- 
ity we used the present performance of one single LIGO in- 
strument and the improved, possible future performance of an 
Advanced LIGO (AdvLIGO) detector (Shoemaker 2007). Note 
that Advanced LIGO possesses several adjustable frequency re- 
sponses. While the 'burst' selection provides a broad range of 
nearly maximum sensitivity (optimal for model sl5g), a 'nsns'- 
tuned instrument is likely to be used for sources that radiate at 
lower frequencies (model sl5h). 

3. Results 

The computed quadrupole wave amplitudes of models sl5g and 
sl5h are displayed in Figs.[3]|4]and|6] A quantitative overview of 
key properties and results is given in Table[T]and Fig.|5] Detector 
dependent quantities are summarised in Table |2] and Fig. [8] 



3.1. Model s15g 

Model sl5g undergoes a rotational core collapse as we assume 
an initial central angular velocity of Q = 27r rad/s. During the 
stage of contraction, the core spins up massively while becoming 
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Table 1. Summary of GW related quantities for the models sl5g and sl5h. yS = is a convenient measure of the average rotation 

rate, which is calculated by taking the ratio of rotational to gravitational energy. Eg„, is the total energy released in gravitational 
radiation and fi, is the frequency of strongest emission at bounce. Finally, we present the maximum amplitudes at different stages of 
their time-evolution in polar (I) and equatorial (II) direction. The subscript / stands for initial, while b and pb stand for bounce and 
post — bounce. 



Model pi 


Pi, 




direction 


|A+.i|[cm] 


|Ax.ft|[cra] 


\A+,pb\[cm\ 


|Ax.pi,|[cm] 


sl5g 0.26 X 10-' 


5.2 X 10^2 


1.52 X 10-^ 


I 
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59 








II 


566 


< 1 


30 


2 


sl5h 0.59 X 10"^ 


1.7 X lO^"* 


4.72 X lO^'" 
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13 
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II 


8 


< 1 
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s15g 
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-0.02 



0.04 



Fig. 2. Time evolution of the central densities of the models sl5g 
and sl5h. At core bounce the following values are assumed: 
Pc - 3.8 X 10'** g/cm-' for model sl5g and pc - 4.5 x 10'"* g/cm-' 
for model sl5h. The density decrease in the postbounce phase is 
most likely due to numerical dissipation of steep density gradi- 
ents in the neutron star. 



oblate. The collapse gets abruptly halted due to the stiffening of 
the EoS above nuclear densities. Afterwards the core rebounds 
and drives a hydrodynamical shock wave outwards. These con- 
ditions give rise to strong time-dependent variations in the phys- 
ical quantities that affect the quadrupole tensor (see equ|7]l. Its 
behaviour is then directly reflected by the model's gravitational 
wave signature. Since the core collapse proceeds nearly axisym- 
metrically, the only GW amplitude considerably driven by the 
rotationally induced large-scale asymmetries is A+//. It exceeds 
the other sl5g-wave trains A+/,Ax/ and Ax// by 1-2 orders of 
magnitude, as one can see in Fig. [3] for times around bounce. 
The initial GW signal is emitted just before core bounce, when 
the rapid infall of matter and the spin-up of the core is dominant. 
It shows a prebounce rise. Then, the instantaneous slowdown of 
matter at core bounce leads to a prominent negative peak, which 
is followed by a ring-down behaviour that lasts for the first few 
ms postbounce. This generic type I wave characteristics is dis- 
played in the lower left panel of Fig. [3] Performing a Fourier 
transform of the GW signal around bounce (-5 ms < f < 5 
ms), we find a spectrum with a very narrow bandwidth peaking 
around 893 Hz. This is in good agreement with the recent find- 
ings from Miiller et al. ( 2004), b ut more than 150 Hz higher than 
in (iDimmelmeier et al.i,2007t lOtt et al.ll2007h . It has been sug- 
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Fig. 3. Time evolution of the quadrupole amplitudes A-hI, AxI, 
A-hII and Axil for model sl5g. 
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Fig. 4. Time evolution of the quadrupole amplitudes A-hI, AxI, 
A-hII and Axil for model sl5h. 



gested that the difiference stems from the fact that full relativistic 
calculations shift the bounce spectrum to lower frequencies in 
comparison to t he one s using an efifective gravitational potential 
(iDimmelmeieil (l2007h V 

After this first and predominantly axisymmetric stage of 
GW-emission, the occurrence of a low r/|W|-instability revives 
the gravitational wave signal again around x 20ms post-bounce 
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dSaiioetalJ (l2003l): IWattsetalJ (l2005h: lOttetalJ (l2005l); 



ISaiio & Yoshida 



(12006^; 'Ou & Tohlin^ (l2006h : lOtt et alJ 
(2007)). Low dynamical instabilities are triggered in 

differentially rotating systems such as neutron stars in situations 
where the patten speed a-p - cr/m of an unstable mode m 
matches the local angular velocity at a point in the star (see 
Fig. Uj, common l y cal led corotation point (the modes are, 
as in IWatts et al.l ( l2005l) . assumed to behave harmonically as 
exp[-!Xcrf - m(^)], where cr is the mode's eigenfrequency). 
It permits the azimuthal fluid modes to amplify. This non- 
axisymmetric process yields a quasi-periodic GW signal with 
a rather constant time-variation, leading to a narrow-band 
emission at 905 Hz which lasts until the end of our simulation, 
as one can see particularly in the upper panels of Fig. [3] for 
times t > 20 ms. The analysis method we use to observe the 
growth of nonaxisymmetric structures decomposes the density 
at a fixed radius R and constant z-component into its azimuthal 
Fourier components as done before e.g. in ref. COu & Tohlin3 
l2006tlOtt et alll2007l) ): 




(16) 



Fig. 5. Time evolution of the GW luminosity in model sl5g. In 
this 70 ms window, about 90% of the energy emitted in GW 
stems from the A-hII contribution around bounce. 



where the complex Fourier amplitudes are defined by 



C, 



1 f^" 
27!- Jo 



p(R,z,<p)e- 



(17) 



In Fig. |6]the normalized mode amplitudes A,„ = \C,„\/Co 
are monitored to measure the growth of unstable modes. In our 
model sl5g we find m = {1, 2, 3)-modes being triggered, with 
the so-called m - 2 bar-mode growing fastest. Further, we state 
that the m - {1,2, 3) modes all possess the same pattern speed. 
The close relation between the m - 2 bar-mode instability and 
the emission of gravitational waves can be seen in the following 
features: First, the sudden onset of GW emission along the pole, 
which must be completely due to nonaxisymmetric dynamics, 
coincides with the amplitude of the m - 2 mode reaching ap- 
proximately the absolute amplitude of the m - A mode caused 
by the grid. Secondly, the dominant frequency of emission cor- 
responds perfectly to the eigenfrequency of the m - 2-mode. 
Finally, the two GW-polarisations + and x are phase shifted 
by ;7r/2, as one would expect of a perfect, monochromatic GW- 
source such as a rotating bar. These findings in the context of the 
low r/|W| instability and supernova d ynamics stand in remark- 
able agreement with the recent ones of lOtt et al] (l2007h . For low 
jS-unstable models similar to sl5g, they found narrow-band GW 
emission at a: 920 - 930 Hz. The main difference to our calcu- 
lations is the point that the dominant mode which was found in 
those computations was the m - 1-mode. As a closing remark to 
model sl5g we state that the time evolution of the energy emit- 
ted by gravitational radiation (see Fig|5]l fits the behaviour of the 
waves. It demonstrates a large peak around bounce at 1.3 x 10^' 
erg/s followed by a ringdown and an oscillating renaissance at 
about 10"^^ - 10"*^ erg/s for times f > 20 ms. 

3.2. Model s15h 

The slow -rotating model sl5h undergoes a quasi-spherically 
symmetric core collapse, consequently showing fairly weak type 
I amplitudes around bounce (see Fig. HJi. On the other hand, 
the marginally present centrifugal forces allow a core bounce 
at higher central densities (4.5 x lO'** g/cm^) than in the previ- 
ously discussed model (3.8x 10'"^ g/cm^) as one can see in Fig.|2] 




X [km] 

Fig. 7. Angular velocity profiles of model sl5g (upper profile) 
and sl5h (lower profile) along the positive x-axis in the equato- 
rial plane at t - ti, x 6ms (S15g) and f = f/, (sl5h), respectively. 
The m = 2 mode pattern speed of sl5g is given by the horizontal 
line. Note first that both models' innermost x 10 km are nearly 
in perfect solid body rotation and secondly that model sl5g is in 
corotation at a radius of about 10 km, spinning with some 2850 
[rad/sec]. 



The general dependence of Ah- II on rotationally induced asym- 
metries gets obvious in the following feature: while the Ah-II 
amplitude from sl5h at bounce is more than one order of mag- 
nitude smaller as in sl5g, all other amplitudes, namely A+I, Axl 
and Axil are of the same order of magnitude as their counterparts 
from model sl5g, emphasising a general feeble subordination to 
rotationally caused aspherities around bounce. 

Prompt convection, based on the negative entropy gradient 
left behind the stalling shock starts just a bit before 10 ms after 
bounce. It leads to a convectively driven rise of the GW sig- 
nal; its characteristic is dominated by the stochastic process of 
convective mass motion: neither a clear signal type nor any cor- 
relation between the two polarisations as can be found here, as 
one expects from this kind of matter motion. However, there is 
a peak in the wave train at a; 7 ms, best visible in A+II, which 
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Fig. 6. The upper panel displays the GW emission of the A+ and the Ax-amplitude along the pole. The lower panel shows the 
normalized mode amplitudes A„, for m - {1,2,3,4) extracted at a radius of 20 km. Note the sudden onset of the nonaxisymmetric 
GW-signal along the pole as soon as the m - 2 mode amplitude reaches approximately the absolute size of the m - 4 grid mode 
background. Notice further that A+ and Ax oscillate at the same frequency of 9Q5Hz, phase shifted by n/l, as one would expect 
from a rotating bar. 



we could not attribute to a physical feature of core collapse. It 
is most likely the result of a grid alignment effect at a radius of 
30-50 km. When the fluid first breaks its spherical symmetry in 
this region, this leads to a numerical quadrupole moment in the 
fluid aligned with the main coordinate axes. The overall spectra 
arising from the sl5h-GW amplitudes is qualitatively and quan- 
titatively rather different to the one of sl5g: it ranges from some 
Hz to about 1000 Hz with major contributions below 500 Hz, 



clearly dominated by post-bounce convective motions (see Fig. 
|8j. The signal lasting from later times (f > 20 ms) mainly takes 
place at a frequency range from roughly 250 Hz down to about 
about 10 Hz. This broad band emission of GW mirrors the wide 
spread and inhomogeneous velocity-distribution of the convec- 
tive motion. The total energy emitted in GW is nearly three or- 
ders of magnitude lower than in the previously discussed model 
(see Table [TJ, and the max luminosity reaches only » 3 x 10"^^ 
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Table 2. Summary of detector-dependent quantities. We assume 
the sources to be located at a distance of 10 kpc as measured by 
one LIGO instrument and at 1 Mpc for Advanced LIGO (sl5g) 
and 100 kpc (sl5h), respectively. The abbreviation li/.??. indicates 
the distance between detector and source, measured in kpc. fc 
and he denote the characteristical frequency and amplitude; SNR 
abbreviates 'signal-to-noise' ratio. We determined the signal-to- 
noise ratios under the assumption of optimal orientation from 
detector and source for polar (I) and equatorial (II) emission di- 
rections, denoted by dir. 



Model 


dist.[k;pc] 


dir. 


MHz] 




SNR 


S^5guG0 


10 


I 


524 


5.1 X lO-^° 


38.4 






II 


469 


6.7 X IQ-^" 


59.1 


slShuGO 


10 


I 


227 


1.4 X 10-21 


3.2 






II 


165 


2.1 X 10-21 


5.6 


S 1 5gAdvLIG0,bursl 


1000 


I 


812 


7.6 X 10-22 


7.3 






II 


841 


1.1 X 10-21 


10.2 


slShAdvLIGO.nsns 


100 


I 


221 


1.5 X 10-22 


3.9 






II 


173 


2.5 X 10-22 


6.6 



erg/s at core bounce. The net size of our convectively driven 
amplitudes (some centimetres) as well as the frequency band of 
emission and the total amount of energy rel eased in gravi t ationa l 
radiation fit qualitatively well the results of iMiiller et al.l (l2004l) . 
where they performed 2D core collapse simulations of slowly 
rotating progenitor models with similar input physics. 

3.3. Detectability 

As in IMiiller et all (|2004 . in Fig. 1 we compare the quantity 
\h(v)\v^^-^ in Eq. dTsl l. evaluated at a Galactic distance of 10 
kpc for a single LIGO detector's /?„,„. The computed SNRs, 
the characteristical frequencies fc and amplitudes he are dis- 
played in Table |2] The results indicate a fair chance of detect- 
ing model sl5g within our Galaxy at the present performance 
of LIGO. Note that the prospect of observing a narrow band 
and long-lasting GW signal from the nonaxisymmetric dynami- 
cal instability (sl5g) is enhanced, increasing the detection-limit 
for a facility running at current sensitivity up to about lOOkpc. 
Furthermore, Fig. |8]reveals that the signal-to-noise ratio of the 
nearly non-rotating supernova model sl5h, where the bounce 
signal is only marginally present, is most probably to small for 
being traced by LIGO, although the peak sensitivity of cur- 
rent detectors lies in the frequency range our results show is 
dominated by convection. The increased detector sensitivities of 
Advanced LIGO should allow the catching of events such as the 
rotational core collapse model sl5g at distances around 1 Mpc 
and the nearly spherically symmetric sl5h at about 100 kpc. 

4. Uncertainties in the underlying models 

Many branches of physics contribute to the models that enable 
the prediction of GW signals from core collapse supernovae. 
Therefore, it is a difficult task to objectively judge on the rela- 
tive importance of uncertainties in GW signal templates. Perhaps 
we can start by distinguishing features that appear in generic 
models from features that appear in deterministic models. A fea- 



ture of generic models would be a feature that appears for all 
representative models out of a given class, while a feature of 
a deterministic model would only appear in a specific simula- 
tion. We believe that most current supernova models should be 
seen as generic models that aim to represent a class of stars, ro- 
tation rates, etc. instead of a specific event. Hence, we should 
concentrate on the generic features of the predicted GW signals 
and defer the discussion of deterministic features to future im- 
proved supernova models. Among the generic features we can 
further distinguish qualitative statements and quantitative state- 
ments. For example, the finding that the A+// quadrupole from a 
rotating collapse model shows a pre-bounce rise before a large 
negati ve peak as shown in Fig. 3 is a generic qualitative state- 
ment (iDimmelmeier et aLll2007h . while the information that the 
Fourier transform peaks at 893 Hz is a quantitative statement. 

It is almost impossible to draw the borderline between 
generic and deterministic features based on a single simulation. 
But even if a full series of models is available, it is difficult to 
exclude that some generic features of this series are not influ- 
enced by a generic limitation of the underlying numerical algo- 
rithms. Hence, it is very important that different groups do not 
just develop the one and best-suited code for a given problem, 
but rather a variety of different numerical approaches so that the 
results can be compared in order to reveal the common generic 
features of the GW signal from core-collapse simulations. 

The last stage of the evolution of a massive star proceeds 
through very diff'erent phases. First, there is the stellar core col- 
lapse and bounce at nuclear density. Then there is a possibly 
extended accretion phase that eventually leads to the supernova 
explosion. These two phases involve different conditions of mat- 
ter and pose different challenges to the numerical modelling. The 
quality and reliability of the models depends very strongly on the 
investigated phase. In the following, we discuss the uncertainties 
of the models for each phase separately. 

4. 1 . Core-collapse and bounce 

The models of stellar core collapse and bounce provide a link be- 
tween progenitor star properties and the first strong emergence 
of a GW signal. Many previous studies analysed the GW signal 
based on idealised input physics in order to study the qualitative 
dependenc e of the core-bounc e signal on progenitor star prop- 
erties (e.g. IZwerger & MiilleJ (119971) ). For example, polytropic 
equations of state are used and the influence of neutrino interac- 
tions on the collapse dynamics are ignored. This approach pro- 
duced a variety of GW bounce signals as a function of param- 
eters like the angular momentum profile of the progenitor star, 
the adiabatic index of the equation of state, or the presence of 
strong magnetic fields. As there is no GW emission in spherical 
symmetry, these studies have always been conducted in a multi- 
dimensional setting. 

However, the input physics relevant during the collapse 
of the stellar core is very rich on the microscopic level 
jMartmez-Pinedo et al.ll2006i) and this dynamical phase is very 
sensitive to small perturbations. The evolution of the collapsing 
core depends significantly on the adiabatic index of the equation 
of state, general relativistic effects, the electron capture rates on 
free protons and nuclei, and the coherent scattering opacities of 
neutrinos off the different nuclei. Input physics improvements 
have been explored in the context of a long-term modelling ef- 
fort to clarify the supernova explosion mechanism in spherically 
symmetric models. Through the ever-increasing power of com- 
puters, these improvements can now be carried over to the multi- 
dimensional predictions of GW signals dMiiller et al.l2004l) . Few 
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Fig. 8. Spectral energy distribution of the GW signals at a distance of 10 kpc in comparison with the present LIGO strain sensitivity 
and the pos sible performance of Advanced LIGO (broadband tuning). The strain sensitivity curves have kindly been provided by 
I Shoemaker! (|2007). The left panel shows the result for the rotational model sl5g and the right panel shows the result for the non- 
rotational model sl5h. Optimal orientation between source and detector is assumed and the polarizations are combined according 
to Eq. ( fTSl ). The dashed lines show the GW signal for an observer on the rotational axis of the source, the solid lines show the 
GW signal for an observer in the equatorial plane. Note model sl5g's signal enhancement of the in the equatorial plane around 
900 Hz finds its physical origin in the core bounce, while the polar signal around the same frequency is dominated by the m = 2 
dynamical instability. Note also the spectral energy distribution of sl5h peaks in the region of maximal detector sensitivity. This 
would probably allow to observe the post-bounce GW-emission from initially nearly spherical symmetric core collapse supernovae. 



groups have implemented equations of state that contain nu- 
clear input physics and added our simple and efficient neutrino 
physics parametrisation scheme to get a more accurate map- 
ping of the progenitor properties onto the expected GW signal 
from the bounce at nuclear matter densities ('Dimm elmeier et al] 
[2007; Cerda-Duran et al. 2007). It was possible to show by two- 
and three-dimensional general relativ istic simulations that only 
type I signals are expected to occur iDimmelmeier et al. I l2007t 
lOtt et alj|2007h . which we independently confirm by the models 
presented in this paper. 

On the other hand, the uncertainties in the progenitor prop- 
erties and the equation of state above nuclear density are rather 
large. Hence, even with the most accurate numerical scheme it 
is not possible to calculate a definitive GW signal from core- 
bounce. The improved models of the latest generation help to 
accurately map progenitor properties and equation of state prop- 
erties to a potentially observable GW signal. This may lead to 
constraints from future GW observations. For example, the am- 
plitude and timing of the GW signal with respect to the neutrino 
signal strongly depends on the rotation rate of the inner stel- 
lar core. It could also depend on the size of the collapsing core 
(weak interactions) and on asymmetric perturbations induced by 



convection in the stellar envelope at the time of collapse. The 
peak frequency of the Fourier transformed signal contains in- 
formation on the compressibility of nuclear matter at bounce. 
Hence, it should be possible to compare GW wave signals for 
different physical equations of state. One should also investi- 
gate a potential impact of strong magnetic fields on the asymme- 
tries of the collapse that might become visible in the GW signal 
jKotake et al.]|2004l;IObergaulinger et alj|2006l) . 

4.2. Postbounce phase 

The homologously collapsed stellar core forms a neutron star 
before it probably enters an extended accretion phase. During 
this postbounce phase layers from the outside of the original 
iron core fall into the standing accretion shock that results from 
core-bounce. They are shock heated and dissociated, and set- 
tle on the protoneutron star. For several reasons, this phase is 
much more difficult to accurately capture in a numerical model 
than the collapse phase. One problem is that the neutrinos do 
not only stream away from the surface of the neutron star, their 
fractional absorption behind the standing accretion shock leads 
to an essential feedback to the accretion dynamics. The hot lay- 



10 



S. Scheidegger et al.: Gravitational waves from 3D MHD core collapse simulations 



ers around the protoneutron star show several three-dimensional 
fluid instabilities that make an accurate treatment of the radiative 
coupling between the deleptonisation of the protoneutron star 
and the accreting layers very delicate. Additionally, one has to 
keep in mind that the density contrast between the center of the 
neutron star and these hot layers behind the standing accretion 
shock may easily exceed five orders of magnitude, which im- 
poses severe time step constraints on codes that treat the whole 
domain consistently for the extended evolution time until the su- 
pernova explosion is thought to be launched. As the time step of 
the hydrodynamics part is limited by the ratio of the zone width 
to the signal speed, larger time steps can be taken if the resolu- 
tion is low in the regions of largest sound speed. As described 
above, we can currently perform simulations with an equidis- 
tant resolution of 1 km. This resolution is conveniently high in 
the outer layers behind the standing accretion shock, but rather 
low at the surface and in the interior of the protoneutron star. 
We avoid numerical artifacts from the steep density gradient at 
the surface of the protoneutron star by analytically considering 
the hydrostatic density gradient in the TVD advection scheme. 
However, the low resolution might still suppress an efficient cou- 
pfing of low protoneutron sta r modes to higher dynamical modes 
jWeinberg & Ouataertll200^ . A more flexible grid for our code 
is under development in order to adapt the resolution to the local 
conditions. 

For these reasons, we consider the technical uncertainties in 
the postbounce phase to be on a similar level than the uncertain- 
ties of the input physics, which of course continue to be present 
after core-bounce. Only very few studies have predicted grav- 
itational waves based on a postbounce model that includes so- 
phistic ated neutrino physics in axisymmetry (e.g. iMiiller et alj 
(12004 ). Neutrino physics approxim ations have been used in 
earlier three-dimensional models e.g. lFrver et alj (|2004|) . In our 
current models, we continue the hydrodynamical simulation to 
the postbounce phase, but are aware that the neutrino physics 
parametrisation cannot handle the neutronisation burst that sets 
in shortly after bounce and it does not feature the deep electron 
fraction trough that develops behind the shock due to the copi- 
ous electron captures. Furthermore, the neutrino parametrisation 
scheme does not feature any neutrino heating. Hence, it is im- 
portant to distinguish the results from the collapse phase and 
bounce, where we believe that the neutrino physics parametri- 
sation is a viable and reasonably accurate approach, from the 
dynamics of the postbounce phase, where the neutrino physics 
parametrisation is not sufficiently accurate. The data from the 
postbounce evolution in our models has to be understood as an 
idealised exploration of potential GW features after the bounce 
signal. This is currently the state-of-the-art in three-dim ensional 
simula tions and handled the same way in the models of lOtt et alj 
(12001) . 

In order to improve the models also in the postbounce 
phase, we have developed the isotropic diff'usi on source approx- 
imati on (IDS A) for the neutrino transport ( Liebendorfer et af] 
l2007h for future models. The scheme has been implemented 
and tested in spherical symmetry, but is not yet fully functional 
in three dimensions. From basic considerations and compari- 
son t o the results in axi symmetry with accurate neutrino trans- 
port (iMiiller et alj|2004i) we expect that the additional neutrino 
emission will lead to a more compact neutron star that is ac- 
creting matter from an envelope with stronger fluid asymme- 
tries. Compared to our current simulations, one could expect a 
stronger GW signal in the late postbounce phase with a shorter 
periodicity in the signal for comparable initial rotation rates of 



the progenitor. However, this will be explored in more detail with 
the next generation of postbounce supernova models. 

5. Conclusion 

In this paper we presented the GW analysis of two 3D MHD core 
collapse supernova simulations, which differ only in the amount 
of initial rotation. Model s 1 5h started with an angular velocity of 
Q. - 0.3 rad/s, model sl5g with Q. - 2n rad/s. We incorporated 
progenitor stars from stellar evolution calculations, spherically 
symmetric GR effects, the Lattimer-Swesty EoS, magnetic fields 
and a neutrino parametrisation scheme that is accurate until the 
first few ms after bounce. The particular choice of the initial 
angular momentum allowed clear distinctions between gravita- 
tional wave features that stem from rotation or nonaxisymmetric 
motions. 

The amplitudes for the three directions and polarisations, 
A+/, Ax/, Ax//, are not very sensitive to rotation and show a sim- 
ilar size of several centimetres. Apparently, they initially cou- 
ple only weakly to rotationally induced large scale asymmetries 
in the mass-energy distribution. The only GW amplitude that is 
strongly correlated to axisymmetric rotation at the time of core- 
bounce turns out to be A+// in the 9 - 7r/2, (p - 0-direction. 
In model sl5g, it exceeds all other amplitudes and in particu- 
lar the corresponding one from model sl5h by more than one 
order of magnitude. It shows a clear type I characteristics and 
implies that a rotational core collapse stays ax isymmetric in th e 
early postbounce phase, as lately discussed in lOtt et akl ( l2007h . 
At the bounce stage of the simulation sl5g, the dominant band 
of emission is peaking around 893 Hz (sl 5g), which is i n goo d 



agreement with the recent findings from IMiiller et al 



but roughly 15 Hz higher than in dDimmelmeier et al 



20041) 
120071: 



lOtt et al.ll2007h . Within the investigated time window of about 
70 ms duration, the channel A+u in model 15g accounts during 
the bounce and ring-down phases for » 90% of the total energy 
release in GW emission. 

The models of the later postbounce phase are yet affected by 
significant technical and physical uncertainties. They should be 
understood as a qualitative outlook to future models that will 
have to include better neutrino transport and a higher resolu- 
tion of the protoneutron star In the postbounce phase, nonax- 
isymmetric dynamics starts to play an important role. The wave 
trains are revived either through convective or through nonax- 
isymmetric instabilities in the protoneutron star. The GW am- 
plitudes from both models show a sustained signal of approx- 
imately constant size, but of different physical origin. In model 
sl5g, the occurrence of a so-called low T l\W\ instability of dom- 
inant m - 2 character around 20 ms post-bounce leads to the 
prolonged narrow band GW emis sion at a frequen cy of 905 Hz, 
which is in good agreement with lOtt et al.l ( l2007b in all points, 
except their dominant mode is m = 1 . On the other hand, the 
onset of nonaxisymmetric dynamics in model sl5h is triggered 
by convective motions due to a negative entropy gradient. In our 
current models, the characteristic frequency of emission linked 
to convective features takes place below 250 Hz. This number 
may require adaption after the inclusion of more accurate neu- 
trino transport. Further, no sign of a decay, be it due to a dynami- 
cal instability or convection, is present in the long term ev olution 
of the w ave trains, a s it was already observed in CMiiller et al.l 
2004; Ot t et al.ll2007l) . 

Finally we conclude that gravitational waves from the dis- 
cussed dynamical features, namely the rotational core bounce 
and its subsequent ringdown and low T l\W\ instability related 
to a core collapse supernova like model sl5g could possibly 



S. Scheidegger et al.: Gravitational waves from 3D MHD core collapse simulations 



11 



be detected by a current LIGO instrument within a distance of 
10 kpc. For model sl5h it seems more likely that LIGO would 
only catch signals from later stages of the supernova evolution, 
meaning frequency distributions caused by convective motion. 
The future improvement of detector sensitivities in combination 
with several adjustable frequency responses of Advanced LIGO 
will give the opportunity to tune the instruments to a particu- 
lar event/source and therefore allow to look much deeper into 
space. This will enhance the chance of observing core collapse 
supernovae similar to our models up to distances of 1 Mpc. 
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